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ABSTRACT 

Context. In gamma-ray binaries such as LS 5039 a large number of electron-positron pairs are created by the annihilation of primary 

very high energy (VHE) gamma-rays with photons from the massive star. The radiation from these particles contributes to the total 

high energy gamma-ray flux and can initiate a cascade, decreasing the effective gamma-ray opacity in the system. 

Aims. The aim of this paper is to model the cascade emission and investigate if it can account for the VHE gamma-ray flux detected 

by HESS from LS 5039 at superior conjunction, where the primary gamma-rays are expected to be fully absorbed. 

Methods. A one-dimensional cascade develops along the line-of-sight if the deflections of pairs induced by the surrounding magnetic 

field can be neglected. A semi-analytical approach can then be adopted, including the effects of the anisotropic seed radiation field 

from the companion star. 

Results. Cascade equations are numerically solved, yielding the density of pairs and photons. In LS 5039, the cascade contribution 
to the total flux is large and anti-correlated with the orbital modulation of the primary VHE gamma-rays. The cascade emission 
dominates close to superior conjunction but is too strong to be compatible with HESS measurements. Positron annihilation does not 
produce detectable 511 keV emission. 

Conclusions. This study provides an upper limit to cascade emission in gamma-ray binaries at orbital phases where absorption is 
strong. The pairs are likely to be deflected or isotropized by the ambient magnetic field, which will reduce the resulting emission seen 
by the observer. Cascade emission remains a viable explanation for the detected gamma-rays at superior conjunction in LS 5039. 

Key words, radiation mechanisms: non-thermal - stars: individual: LS 5039 - gamma rays: theory - X-rays: binaries 



1. Introduction 



The massive star in gamma-ray binaries plays a key role in the 
formation of the very high energy (VHE, >100 GeV) radiation. 
The large seed photon density provided by the O or Be compan- 
ion star, contributes to the production of gamma-rays via inverse 
Compton scattering on ultra-relativistic electrons accelerated in 
the system (e.g. in a pulsar wind or a jet). The same photons 
annihilate with gamma-rays, leading to electron-positron pairs 
production y + y — > <? + + e~. In some tight binaries such as 
LS 5039, this gamma-ray absorption mechanism is strong if the 
VHE emission occurs close to the compact object. Gamma-ray 
absorption can account for the presence of an orbital modula- 
tion in the VHE gamma-ray flux from LS 5039, as observed b y 
HESS dBottcher & Dermerf(2005klBednareld200dlDubusll2006h . 

A copious number of pairs may be produced in the surround- 
ing medium as a by-product of the VHE gamma-ray absorption. 
If the number of pairs created is large enough and if they have 
enough time to radiate VHE photons before escaping, a size- 
able electromagnetic cascade can be initiated. New generations 
of pairs and gamma-rays are produced as long as the secondary 
particles have enough energy to boost stellar photons beyond 
the pair production threshold energy. Because of the anisotropic 
stellar photon field in the system, the inverse Compton radia- 
tion produced in the cascade has a strong angular dependence. 
The cascade contribution depends on the position of the primary 
gamma-ray source with respect to the massive star and a distant 
observer. 



The VHE m odulation in LS 5039 was explained in 
Dub us et all (2008) using phase-dependent absorption and in- 
verse Compton emission, ignoring the effect of pair cascading. 
This model did not predict any flux close to superior conjunc- 
tion, i.e. where the massive star lies between the compact ob- 
j ect and the observer. Th is is contradicted by HESS observations 
dAharonian et alj200 6a). Interestingly, this mismatch intervenes 
at phases where yy-opacity is known to be large r yr » 1 . The 
development of a cascade could contribute to the residual flux 
observed in the system, with secondary gamma-ray emission fill- 
ing in for the highly absorbed primary gamma-r ays. This possi- 
bility w as proposed to explain this discrepancy ( Aharonia n et alj 
2006a) and is quantitatively investigated here in this article. 

The ambient magnetic field strength has a critical impact on 
the development of pair cascading. If the magnetic field strength 
is small enough to neglect the induced deflections on pair tra- 
jectories then the cascade develops along the line of sight join- 
ing the primary source of gamma-rays and a distant observer. 
The particles do not radiate synchrotron radiation. Cascade cal- 
culations are then reduced to a one-dimension problem. Such a 
situation would apply in an unshocked pulsar wind where the 
pairs are cold relative to the magnetic field carried in the wind. 
This paper explores the development of an one-dimensional pair 
cascade in a binary and its implications. 

Previous computations of cascade emis sion in 
binary environment were ca rri ed o ut by Bednarekl 
dl997h: ISierpowska & Bednarekl d2005l): lAharonian et alj 
( 200 6bl): iBednarekl d2006l~ ' 
iKhangulvan et all (|2008); 



2007); Orella na et all (200 
Sierpowska-Bartosik & Torre 
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Fig. 1. This diagram describes the system geometry. A gamma- 
ray photon of energy e\ from the primary source (compact ob- 
ject) interacts with a soft photon of energy eo at a distance r from 
the source and R from the massive star (assumed point-like and 
mono-energetic), producing a pair <? + / e boosted toward a dis- 
tant observer. The system is seen at an angle tfr. 

(2008k IZdziarski etail (120091) . Except for lAharonian etall 
(2006b), all these works are based on Monte Carlo methods. 
One peculiarity of the gamma-ray binary environment is that 
the source of seed photons for pair production and inverse 
Compton emission is the high luminosity companion star. This 
study proposes a semi-analytical model for one-dimensional 
cascades calculations, taking into account the anisotropy in 
the seed photon field. The aim of the paper is to investigate 
and compute the total contribution from pair cascading in the 
system LS 5039, and see if it can account for the measured 
flux close to superior conjunction. The next section presents 
the main assumptions and equations for cascade computations. 
The development and the anisotropic effects of pair cascading 
in compact binaries are investigated. The density of escaping 
pairs and their rate of annihilation are also calculated in this 
part. The cascade contribution along the orbit in LS 5039 is 
computed and compared with the available observations in 
Section 3. The last section concludes on the implications of 
one-dimensional cascades in gamma-ray binaries. More details 
about pair production are available in the appendices. 

2. Anisotropic pair cascading in compact binaries 

2.1. Assumptions 

This part examines one-dimensional cascading in the context of 
binary systems. The massive star sets the seed photon radiation 
field for the cascade. For simplicity, the massive star is assumed 
point-like and mono-energetic. This is a reasonable a pproxima- 
tion as previous studies on absorption (Dubus 2006) and emis- 
sion dDubus et al.l2008i) have shown. The effects of the magnetic 
field and pair annihilation are neglected (see §2.5). Triplet pair 
production (TPP) due to the high energy electrons or positrons 
propagating in a s oft photon field (y + <? + ~ — > e + ~ + e + + e~, 
Mastichia disl 199 lb is not taken into account here. The cross sec- 
tion for this process becomes comparable to inverse Compton 
scattering when E e eo <: 250{m t ,c 2 ) 2 that is for electron ener- 
gies E e ^ 6 TeV interacting with eo ~ 10 eV stellar photons. 
With a scattering rate of about ~ 10~ 2 s _1 , only a few pairs can 
be created via TPP by each VHE electron, before it escapes or 
loses its energy in a Compton scattering. The created pairs have 
much lower energy than the primary electrons. TPP cooling re- 
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Fig. 2. Cascade development along the path to the observer. The 
primary source of photons, situated at the location of the com- 
pact object, has a power law spectral distribution with photon 
index -2 (dotted line). Spectra are computed using the parame- 
ters appropriate for LS 5039 at superior conjunction (d ~ 2/?*, 
R+ = 9.3 R @ , 7V = 39 000 K) for ^ = 30°. The transmitted 
spectrum, including cascade emission, is shown at various dis- 
tances from the primary source: r = R+/4 (black dashed line), 
R-k/2, R+, 2R+ (solid lines) and r = +oo (dotted-dashed line). 
Pure absorbed spectra are shown for comparison (light dashed 
line). 



mains inefficient compared to inverse Compton for VHE elec- 
trons with energy ^ PeV. HESS observations of LS 5039 show 
a break in the spectrum at a few TeV so few electrons are ex- 
pected to interact by TPP in the cascade. Observations of other 
gamma-ray binaries also show steep spectra but this assumption 
will have to be revised if there is significant primary emission 
beyond a; 10 TeV. Pair production due to high energy gamma- 
rays interacting with the surrounding material is also neglected. 
This occurs for y-rays > 1 MeV and the cross-section is of order 
0.04cr T Z 2 cm 2 (see e.g. lLongairlll992l) . with cr T the Thomson 
cross-section. Since the measured Nh is at most 10 22 crrT 2 in 
gamma-ray binaries, pair production on matter will not affect 
the propagation of gamma-rays towards the observer. 

Due to the large velocity of the center-of-mass (CM) frame 
in the observer frame, the direction of propagation of pairs cre- 
ated by yy-absorption is boosted in the direction of the ini- 
tial gamma-ray. For a gamma-ray of energy e\ — 1 TeV, the 
Lorentz factor of the CM to the observer frame transform is 
y' ~ ei/2m e c 2 = 10 6 » 1 (see the appendix, Eq. IA.2I >. Pairs 
produced in the cascade are ultra-relativistic with typical Lorentz 
factor y e ~ 10 6 » 1. Their emission is forward boosted within 
a cone of semi-aperture angle a ~ l/y e <k 1 in the direction of 
electrons. The deviations on the electron trajectory due to scat- 
tering in the Thomson regime are ~ eo/m e c 2 <K l/y e . In the 
Klein-Nishina regime most of the electron energy is given to the 
photon. It is assumed here that electrons and photons produced 
in the cascade remain on the same line, a good approximation 
since y' and y e » 1. This line joins the primary gamma-ray 
source to a distant observer (Fig. U}. 

Pair cascading is one-dimensional as long as magnetic devi- 
ations of pairs trajectories along the Compton interaction length 
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Fig. 3. Spectra as seen by an observer at infinity, taking into account the effect of cascading. Calculations are applied to LS 5039 at 
periastron for different viewing angle if/ = 30°, 60°, 90°, 120° and 150°. Left panel: Complete spectra (solid line) are compared to 
the pure absorbed (light dashed line) and injected (dotted line) spectra. The contribution from the cascade is presented in the right 
panel. 



Ai C remain within the cone of emission of the electrons. This 
condition holds if A ic /(2R L ) < l/y e , with R L the Larmor ra- 
dius. For a typical interaction length Aj C ~ l/(n*cr,- c ) ~ 10 11 cm 
for TeV pairs in LS 5039, the ambient magnetic field must 
be lower than B S= 10~ 8 G. If the magnetic field strength is 
much greater, pairs locally isotropize and radiate in all direc- 
tions. In between, pairs follow the magnetic field lines and 
the dynamics of each pairs m ust be followed as treated in 
Sierpowska & Bednarek ( 2005). The above limit may appear un- 
realistically stringent. However, since deviations and isotropiza- 
tion will dilute the cascade flux, the one-dimensional approach 
can be seen as maximising the cascade emission. More ex- 
actly, this redistribution induced by magnetic deflections would 
decrease the cascade flux at orbital phases where many pairs 
are produced to the benefit of phases where only a few are 
created. Hence, the one-dimensional approach gives an upper 
limit to the cascade contribution at phases where absorption 
is strong. If the flux calculated here using this assumption is 
lower than required by observations then cascading will be un- 
likely to play a role. Finally, one-dimensional cascading should 
hold in the free pulsar win d as long as the pairs move st rictly 
along the magnetic field. InlSierpo wska & Bednarek ( 2005]) and 
ISierpowska-Bartosik & Torres! d2008l) . the cascade radiation is 
comput ed up to the termination s hock using a Monte Carlo ap- 
proach. Sierpowska & Bednarekl J2005) also include a contribu- 
tion from the region beyond the shock. The cascade electrons in 
this region are assumed to follow the magnetic field lines (in con- 
trast with the pulsar wind zone where the propagation is radial). 
There is no reacceleration at the shock and synchrotron losses 
are neglected. In the method expounded here, the cascade radia- 
tion is calculated semi-analytically from a point-like gamma-ray 
source at the compact object location up to infinity, providing the 
maximum possible contribution of the one-dimensional cascade 
in gamma-ray binaries. 



2.2. Cascade equations 

In order to compute the contribution from the cascade, the radia- 
tive transfer equation and the kinetic equation of the pairs have 
to be solved simultaneously. 

The radiative transfer equation for the gamma-ray density 
n y = dN y /dtde\dO. at a distance r from the source is 



dn y 
dr 



dr 



yy 



dr 



dN 
dtde\ 



dE e 



(1) 



where n e = dN e /drdE e dQ. e is the electrons distribution, n* the 
seed photon density from the massive star and dN/dtdei the 
Compton kernel. The kernel is normalised to the soft photon 
density and depends on the energy E e of the electron and the an- 
gle b etween the photon and the direction of motion of the elec- 
tron (Dubus et al. 2008). In the mono-energetic and point-like 
star approximation the stellar photon density can be estimated 
as Li,/4jTcR 2 io, where is the stellar luminosity, eo ~ 2.7 kT+ 
the mean thermal photon energy and R the distance to the mas- 
sive star (see Fig. [TJ. The absorption rate dT yy /dr is given by 
Eq. (IB. 8b . convoluted to the soft photon density. 

The kinetic equation for the pairs is g iven by the following 
integro-differential equation for y e » 1 ( Blumenthal & Gould] 
I^OUZdziarskilligSSHD'Avezac et alj|2007l) 

dn e r Ee , , s , 

— - = -n e (E e ) <P(E e ,E' e )dE' e 

X+OO p 
n e (E' e )<P(E' e ,E e )dE' e +2 J n* g 7y n 7 de u (2) 

where ViEe, E' e ) is the transition rate for an electron of energy E e 
down-scattered at an energy E' e < E e at r. The first two terms on 
the right side of the equation describe the inverse Compton cool- 
ing of pairs, taking into account catastrophic losses in the deep 
Klein-Nishina regime. In this case, most of the electron energy 
is lost in the interaction and the scattered photon carries away 
most of its energy since e\ = E e — E' e * E e . A continuous losses 
equation inadequat ely describes sizeable stochastic losses in a 
single interaction (Blu menthal & Gouldll 1 9701: IZdziarskill 1989t) . 
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Since the inverse Compton kernel gives the probability per 
electron of energy E e to produce a gamma-ray of energy e\ , the 
scattering rate can be rewritten as 



dN 
dtdE' 



(3) 



The expression of dN/dtdE' e is the same as the Compton kernel 
as described before but gives the spectrum of scattered electrons 
instead of the outcoming photon. The first integral in Eq. (f2|l 
is the inverse Compton scattering rate and can be analytically 
expressed as 



*Jm P c- 



9 (E e , E' e ) dE' e = (n c c n*(r) (1 - p e cos O ) , 



(4) 



where /3 e is the electron velocity in the observer frame and cr ic 
is the total i nverse Compton cross-sect ion (for the full expres- 
sion see e.g. [Rybicki & Light marilll979 , Eq. 7.5). The last term 
in the kinetic equation is a source of pairs from yy-absorption 
coupled with the photon density (see the appendices). The pair 
production kernel g yy is normalised to the soft photon density. 

The anisotropic cascade can be computed by inserting 
the anisot ropic kernels for inverse Compton scattering (see 
Eq. A.6 in lDubus etalT 2008) and for pair production obtained 
in Eq. ( IB. 5b in Eqs. |Q]2). The following sections present cas- 
cade calculations applied to compact binaries, using a simple 
Runge-Kutta 4 integration method. It is more convenient to per- 
form integrations over an angular variable rather than r. Here, 
calculations are carried out using i/r r , the angle between the line 
joining the massive star to the observation point and the line of 
sight (see Fig.Q}- 

2.3. Cascade growth along the line of sight 

Figure |2] presents cascade calculations for different distances r 
from the primary gamma-ray source. For illustrative purpose, the 
source is assumed isotropic and point-like, injecting a power- 
law distribution of photons with an index -2 at r — but no 
electrons. The calculations were carried out for a system like 
LS 5039 and for a viewing angle if/ = 30°. In this geomet- 
ric configuration, absorption is known to be strong (r ry » 40 
for 200 GeV photons) and a significant fraction of the total ab- 
sorbed energy is expected to be reprocessed in the cascade, in- 
verse Compton scattering being also very efficient in this config- 
uration. 

Close to the source (r $ d with d the orbital separation), 
absorption produces a sharp and deep dip in the spectrum (light 
dashed line) but the cascade starts to fill the gap (black solid 
line). The angle if/ r increases with the distance r to the primary 
source. Hence, the threshold energy for pair production increases 
as well. Cascading adds more flux to higher energy gamma-rays 
where absorption is maximum. The cascade produces an excess 
of low energy gamma-rays below the minimum threshold energy 
e\ « 30 GeV. Because these new photons do not suffer from 
absorption, they accumulate at lower energies. This is a well- 
known feature of cascading. 



2.4. Anisotropic effects 

This section investigates anisotropic effects in the development 
of the cascade as seen by a distant observer. Cascades are com- 
puted for different viewing angle if/ at infinity, assuming an 
isotropic power-law spectrum for the primary gamma-rays. 



Table 1. Comparison between the radiated power in the cascade 
P r and the absorbed power P a for different viewing angle if/. This 
table provides also the mean energy of pairs (E e ) in the cascade 
at infinity. 
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The left panel in figure [3] shows the complete spectrum tak- 
ing into account cascading (solid line) compared to the pure ab- 
sorbed power-law (dashed line). Due to the angular dependence 
in the pair production process, larger viewing angles shift the 
cascade contribution to higher energies and decrease its ampli- 
tude (Fig. [3] right panel). The cascade flux is small enough to be 
ignored for if/ ^ 150°. 

Three different zones can be distinguished in the cascade 
spectra. First, below the pair production threshold energy, pho- 
tons accumulate in a low energy tail (photon index ss -1 .5) pro- 
duced by inverse Compton cooling of pairs. For if/ ^ 90°, a low 
energy cut-of f is observed due t o the p airs escaping the system 
dBall & Kirkll200d: ICerutti et al.l l2008). This low energy cut-off 
is at about 0.1 GeV for if/ = 30°. The cutoff occurs when the cas- 
cade reaches a distance from the primary source corresponding 
to if/ r » 90°. Then, the electrons cannot cool effectively because 
the inverse Compton interaction angle diminishes and the stel- 
lar photon density decreases as they propagate. For if/ <: 90°, 
particles escape right away from the vicinity of the compan- 
ion star and no tail is produced. Second, above the threshold 
energy, there is a competition between absorption and gamma- 
ray production by reprocessed pairs, particularly for small an- 
gles where both effects are large. Even if cascading increases the 
transparency for gamma-rays, absorption still creates a dip in the 
spectrum. Third, well beyond the threshold energy, absorption 
becomes inefficient. Fewer pairs are created, producing a high 
energy cut-off (« 10 TeV, for if/ = 30°). Klein-Nishina effects 
also contribute to the decrease of the high energy gamma-rays 
production. 

2.5. Escaping pairs 

The spectrum of pairs produced in the cascade as seen at infin- 
ity is shown in figure [4] The density depends strongly on the 
viewing angle as expected, but the mean energy of pairs lies at 
very high energies ((E e ) ^ 100 GeV, see Table [T]). The accu- 
mulation of very high energy particles can be explained by two 
concurrent effects. Far from the massive star (r » d), most of the 
pairs are created at very high energy due to the large threshold 
energy (almost rear-end collision). The second effect is that in- 
verse Compton losses are in deep Klein-Nishina regime for high 
energy electrons. The cooling timescale increases and becomes 
larger than the propagation timescale of electrons close to the 
companion star, producing an accumulation of pairs at very high 
energies. 

The distribution of pairs allows to assess the fraction of the 
total absorbed energy escaping the system in the form of kinetic 
energy in the pairs. This non-radiated power P e can be compared 
to the radiated power released in the cascade P r . Energy conser- 
vation yields the total absorbed power P a = P e + P r . 

The asymptotic radiated power reached by the cascade is 
compared to the total absorbed power integrated over energy in 
Table Q] The fraction of lost energy increases with the view- 
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Fig. 4. Distribution of escaping pairs seen by a distant observer, 
depending on the viewing angle if/ = 30° (dashed line), 60°, 90°, 
120° and 150° (dotted line). The binary parameters are the same 
than in Fig. [3] 



ing angle. In fact, for \p > 90° most of the power remains 
in kinetic energy. Once the electrons are created, only a few 
have time to radiate through inverse Compton interaction. Below 
{ijj < 90°), the radiative power dominates and the cascade is very 
efficient (recycling efficiency up to 80% for if/ - 30°). The cas- 
cade is fully linear, since the power re-radiate d remains muc h 
smaller than the star luminosity P, <k L* dSvenssonl 11987). 
Self-interactions in the cascade are then ne gligible. This i s also 
a consequence of Klein-Nishina cascading (Zdziarski 1988). In 
addition, interactions between particles in the cascade would be 
forcedly rear-end, hence highly inefficient. 

The created positrons will annihilate and form a 5 1 1 keV 
line. However, the expected signal is very wea k. The annihila - 
tion cross-section is cr ~ 0-7-logy/y (see e.g. Longair|[l992f) . 
The escaping positrons have a very high average Lorentz fac- 
tor y k, 10 5 (Tab. 1) so they are unlikely to annihilate within 
the system. They will thermalize and annihilate in the interstel- 
lar medium. Escaping positrons from gamma-ray binaries are 
unlikely to contribute much to the diffuse 511 keV emission. 
The average number of pairs created along the orbit in LS 5039 
(based on the results to be discussed in the following section) is 
N e ~ 5 x 10 35 s _1 . This estimate does not take into account con- 
tributions from triplet pair production or from the pulsar wind 
(for a pulsar injecting pairs with (y e ) ~ 10 s and a luminosity of 
10 36 erg/s, about 10 36 s" 1 pairs are produced). Gamma-ray bi- 
naries have short lifetimes and it is unlikely there is more than a 
few hundred currently active in the Galaxy. Hence, the expected 
contribution is orders-of-magnitude below the positron flux re- 
quired to explain the diffuse 511 keV emission (~ 10 43 s _1 , 
Knodlsed er et al.l2005l) . Even if the positrons thermalize close to 
or within the syste m (because magnetic fie lds contain them, see 
§5) then, following lGuessoum et alj d2006l) . the expected contri- 
bution from a single source at 2 kpc would be at most ~ 10~ 9 ph 
cirT 2 s _1 , which is currently well below detectability. 



Fig. 5. Orbit-averaged spectra in LS 5039 at INFC (0.45 < 
<p < 0.9, grey lines) and SUPC (<£ < 0.45 or ^ > 0.9, black 
lines) an d comparisons with EGRET (dark) and HE SS (light) 
bowties (iHartmanet all 1 19991: lAharonian et al.ll2006ah . Dotted- 
dashed lines represent the primary source of gamma-rays with 
pure absorp tion, injected at r = 0, computed with the model 
described in Dubus et al. (20081) for a mono-energetic and point- 
like star. Dashed lines show the contribution from the cascade 
and thick solid lines the sum of the primary absorbed source and 
the cascade contributions. 
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Fig. 6. Computed light-curves along the orbit in LS 5039, in the 
HESS energy band (flux > 100 GeV). The cascade contribution 
(dashed line) is compared to the primary pure absorbed source 
(dotted-dashed line) and HESS observations. The thick solid line 
shows the sum of both components. 



3. Cascading in LS 5039 

LS 5039 was detected by HESS (lAharonian et al.ll2005l) and the 
orbital m odulation of the TeV gam ma-ray flux was later on re- 
ported in lAharonian et al.l (|2006a). Most of the temporal and 
spectral features can be understood as a result of anisotropic 
gamma-ray absorption and emission from relativistic electrons 
accelerated in the immediate vicinity of the compact object, 



6 



B. Cerutti et al.: One-dimensional pair cascade emission in gamma-ray binaries 



e.g. in the pulsar wind termination shock dDubus et al.l 12008). 
However, this description fails to explain the residual flux ob- 
served close to superior conjunction where a significant ex- 
cess has been detected (b.lcr at phase 0.0+0.05). The primary 
gamma-rays should be completely attenuated. The aim of this 
part is to find if cascading can account for this observed flux. The 
cascade is assumed to develop freely from the primary gamma- 
ray source up to the observer. The contribution of the cascade as 
a function of the orbital phase is also investigated. 

The primary source of gamma-rays now considered is the 
spectrum calculated in [Dubus et alj d2008l) . Figure [5] shows 
phase-averaged spectra along the orbit at INFC (orbital phase 
0.45 < < 0.9) and SUPC (0 < 0.45 or <f> > 0.9) for the 
primary source, the cascade and the sum of both components. 
The orbital paramete rs and the distance (2.5 kpc) are taken from 
ICasares et al.l d2005l) for an inclination i — 60° so t// varies be- 
tween 30° - 150°. The cascade contribution is highly variable 
along the orbit and dominates at SUPC for e\ ~Z 30 GeV, where 
a large pair-production rate is expected. At INFC, cascading is 
small and negligible compared with the primary flux. With pair 
cascading the spectral differences between INFC and SUPC are 
very small at VHE, contrary to what is observed by HESS. In the 
GeV band, cascades contribute to a spectral hardening at SUPC 
close to 10-30 GeV. 

Orbital light-curves in the HESS energy band give a better 
appreciation of the contribution from both components (Fig. [6]). 
The contribution from cascading is anti-correlated with the pri- 
mary absorbed flux. The cascade light-curve is minimum at in- 
ferior conjunction (<p as 0.72). The non trivial double peaked 
structure of the lightcurve at phases 0.85-0.35 is due to compe- 
tition in the cascade between absorption and inverse Compton 
emission. Absorption has a slight edge at superior conjunction 
(0 w 0.06), producing a dip at this phase. Elsewhere, the primary 
contribution dominates over the cascade emission. At lower en- 
ergies (ej < 10 GeV), the cascade contribution is undistinguish- 
able from the primary source. 

In this configuration, the cascade does add VHE gamma-ray 
emission close to superior conjunction but the expected contribu- 
tion overestimates HESS observations. Decreasing the inclina- 
tion of the system does not help: the cascade flux in the TeV en- 
ergy band increases, since the pri mary source is on average more 
absorbed along the orbit (see §3 in Dub usl2006l) . For i ^ 30°, the 
cascade contribution dominates the primary flux at every orbital 
phases in the VHE band. One-dimension cascades can be ruled 
out by the current HESS observations of LS 5039. 

4. Conclusion 

This paper explored the impact of one-dimensional pair cascad- 
ing on the formation of the very high energy radiation from 
gamma-ray binaries in general, LS 5039 specifically. A signif- 
icant fraction of the total absorbed energy can be reprocessed at 
lower energy by the cascade, decreasing the global opacity of 
the primary source. Anisotropic effects also play a major role on 
the cascade radiation spectrum seen by a distant observer. 

A large contribution from cascading is expected in LS 5039, 
large enough that it significantly overestimates the flux ob- 
served by HESS. One-dimensional cascading is too efficient in 
redistributing the absorbed primary flux and can be ruled out. 
However, the fact that it overestimates the observed flux means 
a more general cascade cannot be ruled out (it would have been 
if the HESS flux had been underestimated). If the ambient mag- 
netic field is large enough (B » 10~ 8 G) the pairs will be de- 
flected from the line-of-sight. For B 10~ 3 G the Larmor ra- 



dius of a TeV electron becomes smaller than the LS 5039 orbital 
separation and the pairs will be more and more isotropised lo- 
cally. All of this will tend to dilute cascade emission compared 
to the one-dimensional case, which should therefore be seen 
as an upper limit to the cascade contribution at orbital phases 
where absorption is strong, particularly at superior conjunction. 
Th e initiated cascad e will be three-dimensional as pointed out 
bv lBednarekl dl997l) . Each point in the binary system becomes a 
potential secondary source able to contribute to the total gamma- 
ray flux at every orbital phases. Cascade emission can still be 
sizeable all along the orbit in LS 5039, yet form a more weakly 
modulated background in the light-curve on account of the cas- 
cade radiation redistribution at other phases. The strength and 
structure of the surrounding magneti c field (from both stars) 
has a strong influence on the cas cade dSierp owska & Bednarek 
120051 : iBosch-Ramon et a l. 2008a bj). More realistic pair cascad- 
ing calculations cannot be treated with the semi-analytical ap- 
proach exposed here. Complementary investigations using a 
Monte Carlo approach are needed to better appreciate the cas- 
cade contribution in gamma-ray binaries. 

Finally, the cascade will be quenched if the created pairs lose 
energy to synchrotron rather than inverse Compton scattering. 
This requires ambient magnetic fields B ~Z 5 G, as found by 
equating the radiative timescales for a 1 TeV electron at peri- 
astron in LS 5039. Such ambient magnetic field strengths could 
be reached close to the companion star. In this case an alter- 
native explanation is needed to account for the flux at supe- 
rior conjunction. A natural one to consider is that the primary 
gamma-ray source is more distant to the massive star. The VHE 
source would not be coincident with the compact object loca- 
tion anymore and w ould suffer less from absorption. In the mi- 
croquasar scenario, iBednarekl d2007l) can account for consistent 
flux with HESS observations at superior conjunction if some 
electrons are injected well above the orbital plane (jet altitude 
z > 10 R+). In addition to LS 5039, this possibility was also con- 
sider ed for the system Cyg X— 1 bv lBosch-Ramon et alj (|2008b) 
and lZdziarski et alj d2009l) . 

In practice, reality may consist of a complex three- 
dimensional cascade partly diluted and partly quenched depend- 
ing upon position, angle and magnetic field configuration. 
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production. This work was supported by the European Community via contract 
ERC-StG-200911. 

Appendix A: Pair production 

The main equations for the pair production process are 
bri efly presented her e . Det ai led calculations ca n be f ound 
in iGould & Schrederl d 19671) iBonometto & Reesl d!971l) and 
Bott cher & Schlickeiserl d 19971) . 

A. 1. Kinematics and cross-sections 

The interaction of a gamma-ray photon of energy e\ and a soft 
photon of energy e ( ) in the observer frame leads to the production 
of an electron-positron pair if the total available energy in the 
center-of-mass (CM) frame is greater than the rest mass energy 
of the pair 

2 ei e (1 - cos 6 ) > W e c 4 , (A.l) 

where m e is the electron mass and 6q the angle between the two 
incoming photons in the observer frame. It is useful to define the 
Lorentz invariant s = e\eo (1 - cos #o) /2. Pairs are produced if 
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s > m 2 c 4 and the velocity /? of the electron-positron pair in the 
CM frame is (3 = (1 - m 2 ^ /s) l/2 . 

The differential cross-section dcr-yy / d(J3 cos 6[) in the CM 
frame depends on j3 and the angle ff y between the outcom- 
ing electron-positron pair and the incoming photons. T he ful l 
expression can be found in e.g. Bono metto & Reesl dl97ll) . 
Eq. (2.7). The differential cross-section presents a symmetric 
structure, peaked at cosflj = ±1 and minimum for cos#j = 0. 
Electrons are mostly created in the same and opposite direction 
with respect to the incoming hard photon direction in the CM 
frame. The double peaked structure is enhanced with increasing 
energy (s » m 2 c 4 ) and becomes less pronounced close to the 
threshold (s ~ m 2 c 4 ). The integration over the angles gives the 
total pair production cross-section o- yy , maxim um close to the 
threshold (see Eq. 1 in lGould & Schrederll 19671) . 

The construction of the CM frame with respect to the ob- 
server frame can be simplified if one of the incoming photons 
carries most of the energy. This case is appropriate in the present 
context. For e\ » 6q, the CM frame can be considered as prop- 
agating along the same direction as the high energy photon. The 
velocity of the CM frame in the observer frame can be expressed 
as 



1 - 



4^ 



1/2 



(A.2) 



The total energy of say the electron E e in the observer frame can 
then be formulated using the Lorentz transform from the CM to 
the observer frames 



s l ' 2 + 



i I 1 4\l/2 . 

j8 I s - m e c : I cos ff x 



(A.3) 



providing a relation between E e and cos 9[ . 

A.2. Rate of absorption and pair spectrum kernels 

A gamma-ray photon going through a soft photon gas of density 
dn/dedQ. is absorbed at a rate per unit of path length / 



dr y 



dl 



-JOT 



dn 
dedQ. 



(1 - cos ff) o-yydedQ. 



(A.4) 



The absorption rate gives the probability for a gamma-ray of en- 
ergy €\ to be absorbed but does not give the energy of the pair 

created in the interaction. 

Following IBonometto & Reesl d 197 lb . the probability for a 
gamma-ray of energy e\ to be absorbed between I and / + dl 
yielding an electron of energy between E e and E e + dE e (with a 
positron of energy E e + ~ e\ — E e for e\ » e) is 



8yy 



■ff: 



dn 



(1 - cos 0) ■ 



do- 



rr 



dedQ.. 



(A.5) 



dedQ. dE e 

As w ith anisotropic inverse Compton scattering ( Dubu s~et al.l 
2008), it is useful to consider the case of a monoenergetic beam 
of soft photons. The normalised soft photon density in the ob- 
server frame is 



dn 
dedQ. 



6(e- eo) 6 (cos 6 - cos O ) 5(<p- <po) , 



(A.6) 



where 5 is the Dirac distribution. Injecting Eq. ( IA.6t into 
Eq. (IA.5t gives the anisotropic pair production kernel, a con- 
venient tool for spectral computations. The detailed calculation 
is presented in Appendix B and the complete expression given 
in Eq. ( IB. 5b . The pair production kernel has a strong angular de- 
pendence and a symmetric structure, centered at E e = e\/2 and 



peaked at E e = E ± (see Appendix B, Fig. IB. 11 ). The effect of 
the angle 6q is reduced close to the threshold where the parti- 
cles share equally the energy of the primary gamma-ray photon 
E e ~ E e + ~ £\/2. Far from the threshold, one particle carries 
away almost all the available energy E e ~ e\ . 

The anisotropic kernel integrated over all the pitch angles, 
in the case of a n isotropic gas of photo ns, is consistent with the 
kernel found bv lAharonian et all (11983). Note that a general ex- 
pression for the anisotropic kernel valid beyond the approx ima- 
tion ei » eo is presented in Bottcher & Schlickeiserl (1 1 9971) . 

A.3. Pair density 

The number of pair created per unit of length path and electron 
energy depends on the probability to create a pair and on the 
probability for the incoming gamma-ray to remain unabsorbed 
up to the point of observation so that 



dN e 
dldEe 



[gyy (E e ) + gyy (&1 ~ E e ) } C ^ 



-T„(/) 



(A.7) 



Because of the symmetry in g yy and since electrons and 
positrons cannot be distinguished here, g yy (ei - E e ) = g yy (E e ). 
The integration over electron energy yields 



gyy (E e )dE e ]e- T ^ = 2 



dr- 



dl 



(A.8) 



The total number of pairs produced by a single gamma-ray 
bathed in a soft radiation along the path I up to the distance r 
is then 



Ne(r) = 2(l - £>- T " (r) ) . 



(A.9) 



For low opacity r yy <k 1, pair production is inefficient and the 
number of particles produced tends to « 2r rr . For high opacity 
Tyy » 1, a pair is always created. 



Appendix B: Anisotropic pair production kernel 

This section is dedicated to the calculation of the pair energy 
spectrum produced in the interaction between a single gamma- 
ray photon of energy e\ and a mono-energetic beam of soft pho- 
tons. The general expression in Eq. JA.5t can be reformulated 
using the relativistic invariant s 



8yy 



dn do-y 
e 2 dedQ. dE, 



-dedsd(p. 



(B.l) 



Combining the expression of /3 with the equations Eqs. ( IA.2HA.3l ) 
and defining x = y' 2 , the differential cross-section variables can 
be written as 



/3(x) 



1 



Amzc x 



/3cos8\(x) 



2E e 



1/2- 



The differential cross-section can then be expressed as 

dCTyy dCTyy <1 ^ COS B[) 



(B.2) 



dEe 



d(f3cos0[) dE e 
2 dcr- 



yy 



^(l_l) 1/2 rf(8cosei)' 



(B.3) 
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Fig. B.l. Anisotropic pair production kernel g ry with eo set at 
1 eV for a head-on collision (6\> = n). The kernel is com- 
puted for ei = 265 GeV (dotted line), 300 GeV, 500 GeV, 
1 TeV and 10 TeV (dashed line). The yielding of pairs occurs 
for ei > 260 GeV. 



The complete general formula to compute the spectrum of the 
pair for a non-specified soft radiation field is 



Syy 



4 JIT eV /j _ 



dn 



do- 



ry 



C-9 



A 1 / 2 dedQ. d(J3 cos ff) 



—dedxdfi, (B.4) 



corresponding to Eq. (2.14) in|Bonometto & Reesl d 197 lh . The 
injection of a mono-energetic and unidirectional soft photon 
density (Eq. |A.6t in this last equation yields 



Syy - 



2(1 -juo) dcr. 



\'/ 2 



where fiQ = cos 6q and 
ei 

xo = . 

2e (l - /A)) 



c/(/? cos ff.) 



(j8(xo),j8 cos flj (*>)}, (B.5) 



(B.6) 



This expression is valid for e\ » eo and 5 > m 2 c 4 . The minimum 
and maximum E+ energy reached by the particles is set by 
the kinematics of the reaction and given by 



E 



l ± U-i 

xo 



1/2 



1 - 



4m 2 c 4 x 'l 



1/2 



(B.7) 



Figure IbTI presents the pair production kernel for different in- 
coming gamma-ray energy e\ . 

Note that a kernel can be calculated as well for the absorption 
rate. Injecting Eq. (IA.6t into Eq. (lA.4b is straightforward and 
gives 



<r/r 



yy 



dl 



(1 — COS Of)) <Tyy(J$) 



(B.8) 
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